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ABSTRACT 



Aims. Looking for the primary rainbow in starlight that is reflected by exoplanets appears to be a promising method 
to search for liquid water clouds in exoplanetary atmospheres. Ice water clouds, that consist of water crystals instead 
of water droplets, could potentially mask the rainbow feature in the planetary signal by covering liquid water clouds. 
Here, we investigate the strength of the rainbow feature for exoplanets that have liquid and icy water clouds in their 
atmosphere, and calculate the rainbow feature for a realistic cloud coverage of Earth. 

Methods. We calculate flux and polarization signals of starlight that is reflected by horizontally and vertically inhomo- 
geneous Earth-like exoplanets, covered by patchy clouds consisting of liquid water droplets or water ice crystals. The 
planetary surfaces are black. 

Results. On a planet with a significant coverage of liquid water clouds only, the total flux signal shows a weak rainbow 
feature. Any coverage of the liquid water clouds by ice clouds, however, dampens the rainbow feature in the total flux, 
and thus the discovery of liquid water in the atmosphere. On the other hand, detecting the primary rainbow in the 
polarization signal of exoplanets appears to be a powerful tool for detecting liquid water in exoplanetary atmospheres, 
even when these clouds are partially covered by ice clouds. In particular, liquid water clouds covering as little as 10%- 
20% of the planetary surface, with more than half of these covered by ice clouds, still create a polarized rainbow feature 
in the planetary signal. Indeed, calculations of flux and polarization signals of an exoplanet with a realistic Earth-like 
cloud coverage, show a strong polarized rainbow feature. 
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1. Introduction 

The discovery of the first exoplan et orbiting; a main se- 
quence star almost two decades ago (jMavor fc Quelodll995[ ) 
inaugurated a new era in astronomy. As of today, more than 
700 exoplanets have been detected (source: The extrasolar 
planets encyclopeadia) . Telescope instruments and satellite 
missions, like for exam ple COROT (COny ection, Rotation 
& plane tary Transits') (iBaglin et al.l l20^. NASA's Kepler 
mission (|Koch et al.lll998[) . HA RPS (High Accura cy Radial 
Velocity Planet Searcher) ( e.g. IPepe et al.ll2004[ ). Super- 
WASP (iDeming et al.l l2012D. and, in the ne ar future, GPI 
(Gemini Planet Imager) (jMacintosh et aP l2008fl on the 
Gemini observatory (first on the telescope on the southern 
hemisphere) and SPHERE (Spectro -Polarimetric High- 
Contrast Exoplanet R esearch) (e.g. iDohlen et al.l l2008t 
iRoelfsema et al.l 1201 1[ ) on ESO's Very Large Telescope 
(VLT), to name a few, will rapidly increase the number 
of detected exoplanets. 

The detection methods used and the accuracy of our 
instruments result in most of the exoplanets detected up to 
today being giants, even though in recent years the lower 
mass limit of our detections has been pushed down allow- 
ing for the detection of more than 30 super-Earth planets. 
With an increasing possibility for the detection of the first 
Earth-like planets in the next decade, an important factor 



to consider is how ready our models will be to interpret the 
observations. 

An important factor that needs to be taken into ac- 
count for future efforts to detect signatures of life on other 
planets is the possible inhomogeneity of th e planetary sur- 
face and atmosphere (|Tinetti et al.l l2006l and references 
therein). The existence of continents, oceans and variable 
atmospheric patterns (cloud patches etc), as well as their 
distributions across the planetary surface can have a large 
impact on the observed signal. For this reason, the models 
that we use to interpret the observations should be able to 
handle inhomogeneous planets. 

There exist a number of models th at deal with the 
brightiiess of inhomogeneous exoplanets (iFord et al.l 2001 



Tinettill2006t iMontafies-Rodriguez et al.l 120061: iPalle et al.l 
2008L to name a few). All of the models show the im- 



portance of planetary inhomogeneity and temporal vari- 
ability o n the mode l led p lanetary signal. The pioneering 
work of iFord et aD (|200ll ) as well as later studies (e.g. 
lOaklev fc Cashll2009f ). show a clear diurnal variability of 
the modelled Earth-as-an-exoplanet signal due to areas of 
different albedo passing in and out of the observational field 
of view. 

Among the factors that influence the planetary signal, 
clouds have a prominent role. Observations of earthshine 
for example, have shown that clouds can induce a consid- 
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erable daily variation in the planetary signal. The amount 
of variability observed differs slightly among the observa- 
tions (- 10% for Palle et al. (200 4), ~ 5% for.Goode e t aL, 
(|200lD and a few percent for iCowan et al.l (|2009D l Cloud 
coverage and variability can also influenc e to a large de- 
;ree t he interpretation of the observations. lOaklev fc CashI 
2009t) for example find that mapping the planetary sur- 
face is only possible for cloud coverages smaller than the 
mean Earth one. Even in the case of giant planets or dwarf 
stars, clouds play a crucial role in defi ning the atmospheri c 
thermal profile and eventually spectra (jMarlev et al.ll2010h . 

At present, the characterization of exoplanets is 
mainly done using planetary transits with instruments 
on e.g. t he Hubble and Spit ze r Space Telescopes 



see e.g. |Ehrenreich et al.l 120071: iDeming et"all 120111: 



iTinetti fc GrifEthll2010t) 7 and for some planets, eveii with 



ground-based ins t rumentation (see e.g. [Snellen et al]|2010t 
iBrogi et al.l l2012t Ide Mooii etal] 12012'). With the transit 
method though, during the primary transit, the observed 
starlight has only penetrated the upper layers of the plan- 
etary atmosphere. Cloud layers at lower altitudes could for 
example block out the signal from lower atmospheric lay- 
ers or a possible planetary surface. Even with the help of 
secondary transits, the characterisation of Earth-like ex- 
oplanets in the habitable zone of a solar-type star would 
not be possible, since these planets would yield too weak a 
signal (Kaltenegger & Traub 2009). Direct observations of 
reflected starlight from the planet could solve this problem, 
since then information from the lower atmospheric layers 
and surface could survive in the observed planetary signal. 
The combination in these cases of flux and polarization ob- 
servations could provide us with a crucial tool to break 
any possible retrieval degeneracies (for example, such as 
between optical thicknesses and single scattering albedo's 
or cloud particle sizes) that flux only measurements may 
present. A first dete ction of an exo p lanet using polarimetry 
was claimed bv iBerdvugina et"aI1 (12008). Su bsequent po- 
larization observations by Wiktorowicd ( 20091) could, how- 
ever, n ot confirm this detection. In 2011 Berdvugina et al] 
(|201lt) presented another detection of the same planet at 
shorter wavelengths, which still awaits confirmation by 
follow -up observations . Telescope instruments like GPI 



(iMacmtosh et al.l I2008D an d SPHERE (e.g. iDohlen et all 
l2008t Roelfsema et al.ll2011l) , which have polarimetric arms 
that have been optimized for exoplanet detection, are ex- 
pected to detect and characterize exoplanets with polarime- 
try in the near future. 

The power of polarization in studying planetary atmo- 
spheres and surfaces has been shown multiple times in the 
past through observations of Solar System planets (includ- 
ing Earth itself) (see for example Ha nsen & Hovenier 1974; 
Hansen fc Travis! 11974 iMishchenkoi Il990i: iTomasko et al.l 



20091) ■ as well as by modeling of sol ar sy stem plane ts or gi- 
ant and Earth-like exoplanets (e. g. |Stam 2003; Sta ni_et aD 



2004; 'Saar fc Seagei^ 1 20031: ISeager et al.li2000c iStamll2008t 
Karalidi et al. 2011). Polarization provides us with a unique 
tool for the detection of liquid water on a plan etary at- 
mosphere and surface. IWilliams fc Gaidoi ([200^ e.g. use 
polarization to detect the glint of starl ight reflected on liq- 
uid s urfaces (oceans) of exoplanets and lZugger et al.l (|2010l 
1201 it) conclude that the existence of an atmosphere as thick 
as Earth's would hide any polarization signature of the un- 
derlying oceans. Probably the most interesting feature to 



look for in the polarization signal of exoplanets is the rain- 
bow. 

The rainbow is a direct indication of the presence 
of li quid water d r oplets in a planetar y atmosphere (see 
e.g. iBailevI l2007t iKaralidi et al] 1201 It) . Its angular po- 
sition depends strongly on the refractive index of the 
scattering particles and slightly on their effective radius 
(see .Karalidi et al., 2011,1 , and references therein) . Its exis- 
tence, most pronounced in polarization observations , can b e 
masked by the existence of ice clouds (jGoloub et al.ll2000[ ). 
The latter are often located above thick liquid water clouds 
and can dominate the appearance of the reflected polariza- 
tion signal for ice c loud optical thicknesses larger than 2 
(iGoloub et al.ll2000l) 

In IKaralidi et al.l ()2011l) . we calculated flux and polar- 
ization signals for horizontally homogeneous model plan- 
ets that were covered by liquid water clouds. The polar- 
ization signals of these planets clearly contained the signa- 
ture of the rainbow. The polarization signals of the quasi 
horizontally inhomogeneous planets (where weighted sums 
of horizontally homogeneous planets are used to approx- 
imate the signal of a horizon tally inhomogeneous planet) 
as presented by IStamI ()2008[ ) also show the signature of 
the rainbow. However, to confirm that realistically hor- 
izontally inhomogeneous planets, with patchy liquid wa- 
ter clouds, and with patchy liquid water clouds covered 
by patchy ice clouds, also show the signature of the rain- 
bow, the algori thm for horizontal l y inhomo geneous planets 
as described in IKaralidi fc Stainl ( accepted! ) is required. In 
this pa per we wi l l do exactly that: using the algorithm of 
Karali di fc Stand (jaccept ed) to look for the rainbow on hor- 
izontally inhomogeneous planets, that are covered by dif- 
ferent amounts of liquid and icy water clouds. 

This paper is organized as follows. In Sect. |2l we give a 
short description of polarized light and our radiative trans- 
fer algorithm, and present the model planets and the clouds 
we use. In Sect. |3l we investigate the infiuence of the liquid 
water cloud coverage on the strength of the rainbow fea- 
ture of a planet in flux and polarization. In Sects. 01 and |S] 
we investigate the infiuence of different cloud layers, re- 
spectively with different droplet sizes and with different 
thermodynamic phases (liquid or ice), on the strength of 
the rainbow feature. An interesting test case for the detec- 
tion of a rainbow feature is of course the Earth itself. In 
Sect, ini we use realistic cloud coverage, optical thickness 
and thermodynamic phase data from the MODIS satellite 
instrument to investigate whether the rainbow feature of 
water clouds would appear in the disk integrated sunlight 
that is refiected by the Earth. Finally, in Sect. [71 we present 
a summary of our results and our conclusions. 



2. Calculating flux and polarization signals 

2.1. Defining flux and polarization 

We describe starlight that is refiected by a planet by a flux 
vector -kF, as follows 



TtF = TT 



F 

Q 

u 

V 



(1) 



where parameter ttF is the total flux, parameters nQ 
and ttU describe the linearly polarized flux and parameter 
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ttV the circularly polarized flux (see e. g. iHansen fc Travis! 
11974 iHovenier et all 120041: IStaini l2008ll . All four parame- 
ters depend on the wavelength A, and their dimensions are 
W m~^m~^. Parameters wQ and irU are defined with re- 
spect to a reference plane, and as such we chose here the 
planetary scattering plane, i.e. the plane through the center 
of the star, the pla net and the observer. Parameter ttV is 
usually small ([Hansen fc Travi s 1974), and we will ignore 
it in our numerical simulations. Ignoring ttV will not intro- 
duce si gnificant errors in our calculated total and polarized 
fluxes ijstam fc Hovenieij r2005l. 

The linearly polarized flux, irFp of flux vector ttF is 
independent of the choice of the reference plane and is given 

by 

ttFp = n^Q^TlP, (2) 

while the degree of (linear) polarization P is defined as 
the ratio of the linearly polarized flux to the total flux, as 
follows 



F F 



C/2 



(3) 



For a planet that is mirror-symmetric with respect to the 
planetary scattering plane, parameter U will be zero. In this 
case we use the signed degree of linear polarization Pg = 
—Q/F, which includes the direction of polarization: if Pg > 
(Pg < 0), the light is polarized perpendicular (parallel) 
to the plane containing the incident and scattered beams 
of light. A planet with patchy clouds will usually not be 
mirror-symmetric with respect to the planetary scattering 
plane, and hence U will usually not equal zero. 

The flux vector ttF of stellar light that has been re- 
flected by a spherical planet with radi us r at a distance d 
from the observer [d ^ r) is given by (|Stam et al.ll2006D 



7rF(A,a) = i^5(A,a)7rFo(A). 



(4) 



Here, A is the wavelength of the light and a the planetary 
phase angle, i. e. the angle between the star and the ob- 
server as seen from the center of the pla net. Furthermore , 
S is the 4x4 planetary scattering matrix (|Stam et al.ll2006[ ) 
with elements Uij and ttFq is the flux vector of the incident 
stellar light. For a solar type star, the stellar flux can be 
considered to be unpolarize d when integrated over the stel- 
lar disk (Kemp et al.lll987D . 

We assume that the ratio of the planetary radius r and 
the distance to the observer d is equal to one, and that the 
incident stellar flux ttFq is equal to 1 W m^^. The 
hence normalised flux ttF^ that is reflected by a planet is 
thus given by 



7rPn(A, a) = -aii(A,a) 



(5) 



(see IStariil 120081 : iKaralidi et al.|[20TTl ). and corresponds to 
the planet's geometric albedo Aq when a = 0°. The corre- 
sponding normalized polarized flux is given by 



1 



7rPn,p(A,a) = -v/a2i(A,a)2 +a3i{X,a)^, 



(6) 



with a2i and asi elements of the planetary scattering ma- 
trix. Our normalized fluxes ttF^ can straightforwardly be 
scaled for any given planetary system using Eq. S] and in- 
serting the appropriate values for r, d and ttPq- The degree 
of polarization P is independent of r, d and ttFq, and will 
thus not require any scaling. 



2.2. The radiative transfer calculations 

The code that we use to calculate the total and polar- 
ized fluxes of starlight that is reflected by a model planet 
fully includes single and multiple scattering, and polariza- 
tion. It i s based on the same efficient adding-doub l ing a l- 
gorithi n Ide Haan et al.l ll987D used bv IStam et al] ()2006[ ): 
IStaml (|2008D . To calculate flux and polarization signals 
of horizontally inhomogeneous planets, we divide a model 
planet in pixels with a size such that we can assume that 
the surface and atmospheric layers are locally plane paral- 
lel and horizontally homogeneous. We then use the code 
for horizontally i nhomogen eous planets, as presented in 
Karalidi fc StanJ ([accepted'): the contribution of every illu- 
minated pixel that is visible to the observer to the planet's 
total and polarized flux are calculated separately, the polar- 
ized fluxes are rotated to the common reference plane, and 
then all fluxes are summed up to get the disk-integrated 
planetary total and polarized fluxes. From these fluxes, the 
disk-integrated degree of linear polarization P is derived 

(Eq.[3l). 

As before ([Karalidi fc StamI lacceotedh . we divide our 
model planets in plane-parallel, horizontally homogeneous 
(but vertically inhomogeneous) pixels of 2° x 2°. Our pixels 
are large enough to be able to ignore adjacency effects, i.e. 
light that is scattered and/or reflected within more than 
one pixel (e.g. light that is reflected by clouds in one pixel 
towards the surface of another pixel). These effects, that 
make the fluxes emerging from a given type of pixel depen- 
dent on the properties of the surrounding pixels, show up 
for higher spatial resolutions, for example, with pixels that 
are smaller than about 1x1 km^ (jMarshak et al.ll2008l ). 



2.3. The model planets 

All model planets have a vertical inhomogeneous atmo- 
sphere on top of a black surface. The assumption of a black 
surface is a very good approximation for an ocean sur- 
face (without the glint), which covers most of the Earth's 
surface. We'll discuss the effects of brighter surfaces on 
our results where necessary. All model atmospheres have a 
pressure and temper ature profile representati ve for a mid- 
latitude atmosphere (jMcClatchev et al.l[l972h (see Fig. [T|). 
We divide each model atmosphere in the same 16 layers. 
The total gaseous (Rayleigh) scattering optical thickness 
of the atmosphere is 0.097 at A = 0.550 fj,m and 0.016 at 
A = 0.865 /.tm. We do not include gaseous absorption, which 
is a good assumption at the wavelengths of our interest. 

The flux and polarization signals of our model planets 
are fairly insensitive to the pressure and temperature pro- 
files, but they are sensitive to the horizontal and vertical 
distribution of the clouds, and to the microphysical proper- 
ties of the cloud particles. We will use two types of clouds: 
the first type consisting of liquid water droplets, and the 
second type consisting of water ice particles. The details 
of these particles will be discussed in Sect. 12.41 The liquid 
water clouds will be located below an altitude of 4 km (i.e. 
at pressures lower than 0.628 bars, which corresponds to 
temperatures higher than 273 K) and the ice clouds above 
that altitude. 

We create our cloud maps by using an ISCCP yearly 
cloud map which we filter according to the clouds' opti- 
cal thicknesses. Figure [5] shows a sample map for a model 
planet with a single cloud layer that covers ~25% of the 
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Fig. 1. The atmospheric pressure and altitude as functions 
of the temperature for the model planets. The liquid water 
clouds are located below 4 km (for example in between 3 km 
and 4 km as in the gray, dotted area), and the ice clouds 
between 7 km and 8 km (blue, dashed-dotted area). 




Fig. 2. A sample map of our model planets. Here, the plan- 
etary atmosphere contains a single layer with clouds that 
cover ~ 25% of the planet (the dark regions are the clouds). 



planet. As shown in iKaralidi fc Staml (|accepted| ) . the pre- 
cise location of clouds can influence in particular the polar- 
ized phase function of a planet (i.e. the degree of polariza- 
tion of the reflected starlight as a function of the planetary 
phase angle). To avoid changes in Tri^n p and P due to the 
location of clouds when the coverage increases, we increase 
the cloud coverage of a planet by letting existing clouds 
grow in size. Cloudy regions thus remain cloudy, while sur- 
face patches are getting cloudy. 

2.4. The cloud particles 

The liquid water cloud particles are spherical, with the 
standard size distribution given by (|Hansen fc Travislll974[ ) 



l(r) — (77-(1^3uoff)/fctf g-r/roftUotf 



(7) 



with n{r)dr the number of particles per unit volume with 
radius between r and r + dr, C a constant of normalization. 



Tcff the effective radius and fcff the effective variance of the 
distribution. Terrestrial liquid water clouds are in general 
composed of droplets with radii ranging from ^ 5 /im to 
~ 30 ^m (|Han et allfTool . We use either small droplets 
(type A), with reS = 2.0 /j,m and v^s = 0.1, or larger 
drop lets (type B) with r ^s = 6.0 /im and Vcs = 0.4 (see 
also 'Karalidi et al."201l'). The latter are similar to those 
used by van Dicdcnhovcn et al.l ()2007t ) as average values for 
terrestrial water clouds. We use a wav elength independent 
refractive index of 1.335-FO.OOOOli fsee lKarahdi et al.ll201lL 
and references therein). For a given wavelength, and values 
of Toff and Woff, we calculate the extinction cross-section, 
single scattering albedo and the single scatterin g matrix 
using Mie-theory (|de Rooii fc van der Stad 1984 ) and n or- 
malized according to Eq. 2.5 of Hansen fc Travis (|1974f) . 

Ice crystals in nature present a large variety of shapes, 
depending on the t emperature and humidity conditions 
during their gro wth (iHessI 19981) . Even though the variety 
of shapes is large. [Maeono fc Lea (|1966[ ) showed that only a 
small number of classes suffice for the categorization of most 
natural ice crystals. Until recently most researchers used 
perfect hexagonal columns and plates in order to model the 
light scattered by natural ice crystals, including halo phe- 
nomena. However, when these particles are o riented, they 
give rise to halos, which are rarely observed (jMacke et al.l 
[1996 : Hess 1998) . For this reason a number of models were 
created to i nodel the signal of ice crystals without halo 
phen omena (|Hess fc Wiegnerj|1994HMacke et aIlll996l:lHessl 
I1998D . In t his pap er we use an updated version of the model 
crystals of iHessI ('1998V Finally, we should mention that 
the sizes of most atmospheric ice crystals are considerably 
larg er than the visible wavelengths we will use in this pa- 
per (jMacke et al.lll996h . In particular, the crystals we use 
in this paper range in size fr om ^6 jum up to ^2 mm, with 
a size distribution based on ([Hevmsfield fc Plattj|1984l ). 

Figure [3] shows the total flux, the polarized flux and the 
degree of linear polarization of light that has been singly 
scattered by the cloud droplets and the ice particles at a 
wavelength A equal to 0.550 /im, when the incident light 
is unpolarized. The curves for scattering by gas molecules 
(Rayleigh scattering) have also been added. All curves have 
been plotted as functions of the phase angle a instead of 
of the more common single scattering angle O to facilitate 
the comparison with the signals of the planets. Since for 
single scattering by the spherical cloud particles or by the 
ensemble of randomly oriented ice crystals, the scattered 
light is oriented parallel or perpendicular to the scattering 
plane (the plane containing the directions of propagation 
of the incident and the scattered light), we use the signed 
degree of linear polarization Ps in Fig. [31 when Pg is larger 
(smaller) than 0, the direction of polarization is perpendic- 
ular (parallel) to the scattering plane. The absolute value 
of Ps equals P. 

The curves for Rayleigh scattering by gas clearly show 
the nearly isotropic scattering of the total flux, the sym- 
metry of the polarization phase function and the high po- 
larization values around a = 90°. The curves for the to- 
tal flux scattered by both types of liquid water droplets 
have a strong forward scattering peak at the largest val- 
ues of a (the smallest single scattering angles), which 
is due to refraction and depends mostly on the size of 
the scattering particles, and not so much on their shape 
(jMishchenko et aT]|2010l ). Another characteristic of the flux 
phase function of the droplets is the primary rainbow, that 
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Fig. 3. Single scattering phase function F, hnearly polar- 
ized flux Q, and degree of linear polarization Pg as func- 
tions of the phase angle a {a = 180° — 8, with Q the 
single scattering angle). Linearly polarized flux U equals 
zero, and Pg is defined as —Q/F. Curves are shown for 
liquid water droplets type A {r^s = 0.2 /im, VeS = 0.1), 
type B (reff = 6.0 nm., = 0-4), an d water ice parti- 
cles (|Hess fc Wieeneil 11994 iHessI I1998D . all calculated at 
A = 0.550 iim. For comparison, the curves for Rayleigh 
scattered light are also shown. 

is due to light that has been reflected once inside the par- 
ticles. For a given wavelength, the precise location of the 



rainbow depends on the particle size: for the small type A 
particles, it is found close to a = 30°, while for the larger 
type B particles, it is close to a = 40°. The primary rain- 
bow also shows up in the scattered polarized flux and very 
strongly in the degree of polarization Pg- The direction of 
polarization across the primary rainbow is perpendicular 
to the planetary scattering plane. The polarized flux and 
Ps go through zero, thus change direction, a few times be- 
tween phase angles of 0° and 180°. These particular phase 
angles are usually referred to as neutral points and, like the 
rainbow, they depend on the particle properties and the 
wavelength. At 0.550 /im, the type A particles have neutral 
points at 5°, 20°, 76° and 158°, and the type B particles at 
2°, 22°, 94°, and 160°. 

The total flpc s c attered by the ice particles 
(jHess fc Wiegneil 119941: iHessI 119981 ) has a smooth ap- 
pearance without a feature such as the rainbow but with a 
strong forward scattering peak at the largest values of a. 
The polarized flux and the degree of polarization of the ice 
particles are also smooth functions of a. The neutral point 
of the ice particles is around a = 140°. 

Fig. |4] shows the degree of linear polarization Pg for 
the type B liquid water droplets and the ice particles at 
A = 0.550, 0.660 and 0.865 /zm. The ice particles are rel- 
atively large, and therefore their scattering properties are 
virtually insensitive to A across the wavelength region of our 
interest (the visible). The scattering properties of the liq- 
uid cloud particles vary slightly with the wavelength: with 
increasing A, the strengths of the primary and secondary 
rainbows decreases slightly, and the neutral point around 
intermediate phase angles shifts towards smaller phase an- 
gles. With increasing A, the primary rainbow also shifts to- 
wards smaller phase angl es (larger single sc attering angles) . 
As also discussed in Kar alidi et al.l (|201lt) . the dispersion 
of the primary rainbow depends on the size of the particles 
that scatter the incident starlight: in large particles, such 
as rain droplets, the dispersion shows the opposite behav- 
ior, with the primary rainbow shifting towards larger phase 
angles with increasing A. This latter shift gives rise to the 
well-known primary rainbow seen in the rainy sky with the 
red bow on top and the violet bow at the bottom. 

3. The influence of the liquid water cloud coverage 

In this section, we explore the strength of the primary rain- 
bow feature as a function of the liquid water cloud cover- 
age. All clouds have the same optical thickness, i.e. 2.0 (at 
0.550 /im), and the same altitudes of their bottoms and 
tops, i.e. 3 and 4 km, respectively (see Fig. [1]). The clouds 
consist of type A droplets. 

Figure [5] shows the total reflected flux irFn, the polar- 
ized reflected flux TrFn.p , and the degree of polarization P 
at 0.550 /xm, as functions of the phase angle a and the per- 
centage of cloud coverage. Horizontal cuts through Fig. [S] 
at a number of coverages are shown in Fig. [S] Because the 
model planets are not mirror-symmetric with respect to 
the planetary scattering plane, we use Eq. |3] to deflne the 
degree of polarization instead of the Pg used in the previous 
section. 

Clearly, as the cloud coverage increases, the flux and 
polarization features converge smoothly towards the signals 
for completely cloudy planets (see also KaialidLgt.al. 20lTl 
for flux and polarization signals of completely cloudy plan- 
ets). At small coverages, P has a strong maximum around 
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Fig. 4. Singly scattered F, Q, and Pg as functions of a 
for the liquid water droplets type B and the ice particles 
at A = 0.550, 0.660 and 0.865 /.tm. The curves for the ice 
particles overlap each other. 
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Fig. 5. The total flux 7rF„, the polarized flux nF^^p, and 
the degree of polarization P as functions of the phase angle 
a for liquid water cloud coverages from 0% up to 100%. 
Wavelength A is 0.550 /im, the optical thickness of the 
clouds is 2.0, and their top is at 4 km. The clouds are com- 
posed of type A droplets. 



a = 90°, which is due to scattering by the gas molecules 
above and between the patchy clouds. The Rayleigh scat- 
tering optical thickness above the clouds is 0.06 at A = 
0.550 /J,m, and between the clouds, it is 0.097 (see Sect. l2.3|l . 
As expected from the single scattering curves (Fig. [S]), the 
primary rainbow feature is located close to a = 30°. The 
precise rainbow phase angle depends on the size of the par- 



ticles and the wavelength, but, as expected, it does not 
depend on the cloud coverage. 

In the total flux TrPn, the rainbow is difficult to discern 
regardless of the cloud coverage (Fig. [SJ, as can also be 
seen in Fig. [51 In the polarized flux, the rainbow shows up 
as a local maximum for coverages of about 30% or more. 
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Fig. 6. Horizontal cuts through the panels of Fig. [5] for 
cloud coverages of 20%, 30%, 40%, and 60%. 



In the degree of polarization, P, the rainbow feature is a 
shoulder on the Rayleigh scattering maximum from a cloud 
coverage of ~20%. For a cloud coverage larger than about 
40%, the rainbow causes a local maximum in P, because 
the Rayleigh scattering maximum decreases with increasing 
cloud coverage. A reflecting (i.e. non-black) surface under- 
neath our atmosphere would increase the contrast of the 
(primary) rainbow peak by lowering the intensity of the 
Rayleigh scattering peak due to the increase of light with 
a generally low degree of polarization. 



At phase angles near 20°, nFa^p and P are close to zero, 
fairly independent of the cloud coverage (see Fig. [6]). This 
phase angle corresponds to a neutral point in the single 
scattering polarization phase function (Fig. [3]). The type A 
cloud particle's neutral point near 76° is lost in the Rayleigh 
scattering contribution to P. The neutral point near 5° 
yields the near-zero polarization region at the smallest 
phase angles, while the neutral point near 158° and the 
generally low polarization values at those phase angles (see 
Fig. [3]) show up as the broad low P region at the largest 
phase angles in Fig. [S] Note that at these largest phase an- 
gles, the planet is almost in front of its star, and will be 
extremely difficult to detect directly anyway. 

Figure [7] shows the same as Fig. [SJ except for A = 
0.865 /im. At this wavelength, the Rayleigh scattering op- 
tical thickness above the clouds is 0.01, while between the 
clouds it is 0.016 (see Sect. 12. 3p . The optical thickness of 
the clouds at 0.865 /im is 2.1. The differences in reflected 
fluxes and the degree of polarization between Fig. [7] and 
Fig. [S] are due to the difference in Rayleigh optical thick- 
ness above and between the clouds and to the difference in 
scattering properties of the cloud particles at the different 
wavelengths. Note that the effects of lowering the Rayleigh 
scattering optical thickness above the clouds are similar to 
increasing the altitude of the clouds. 

At 0.865 /im, the reflected total flux irFn is a very 
smooth function of a, almost without even a hint of the 
rainbow. This smoothness is due to the smoother single 
scattering phase function of the small type A particles at 
this relatively long wavelength. The degree of polarization 
shows a very strong rainbow signature, especially for cloud 
coverages larger than 10%. The reason that the rainbow 
is so strong is mostly due to the small Rayleigh scattering 
optical thickness above the clouds, which suppresses the 
strong polarization maximum around 90°. The strength of 
the rainbow appears to be fairly independent of the percent- 
age of coverage (all clouds have the same optical thickness) . 
Around a = 108°, P is virtually zero for a coverage of 40%. 
These low values of P (also in the graphs for the 20% and 
30% coverage, but then around a = 135° and 130°, respec- 
tively) are due to the direction of polarization of the light 
scattered by the cloud particles (cf. Figs. [3] and H]) that is 
opposite to that of light scattered by the gas molecules. 
Because at 0.865 /im, the Rayleigh scattering optical thick- 
ness is small, the single scattering polarization signatures 
of the ice particles dominate the polarization signature of 
the planet. 

4. The influence of mixed cloud droplet sizes 

In the previous section, all clouds were composed of the 
same type of liquid water cloud droplets. In reality, cloud 
droplet sizes vary across the Earth. Typically, droplets 
above continents are smaller than those above oceans, 
due to different types of condensation nuclei and differ- 
ent amounts of condensation nuclei. In particular, typical 
cloud particle radii range from about 5 /zm to 30 ^m, with 
a global mean value of about 8.5 /im abo ve continental ar - 
eas, and 11.8 /tm above maritime areas (iHan et aLlliOOl . 
Also within clouds, droplet size distributions will vary, de- 
pending e.g. on internal convective updrafts or downdrafts 
that influence dr oplet growth through condensation and 
collisions (see e.g. [Stephens fc PlattlllQSfl : ISpinhirne et al.l 
[19891). Here, we will present the influence of clouds with 
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Fig. 7. Similar to Fig. El except for A = 0.865 /xm. 

different liquid water droplet size distributions on the flux 
and polarization signals of a model planet. 

Our model planet has two layers of clouds, with the 
lower clouds located between 1 and 3 km (0.902 and 
0.710 bar) and the upper clouds between 3 and 4 km 
(0.710 and 0.628 bar) (cf. Fig.[T]). The lower clouds are com- 
posed of the type B droplets (rcs = 6.0^m, VcS = 0.4) and 
have a coverage of 42.3%. The upper clouds are composed 
of the smaller, type A droplets {rcs — 2.0^m, VcS = 0.1) 
and have a coverage of 16%. The maps for the upper and 
the lower clouds were generated separately and hence there 
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Fig. 8. irFn, ttF^.p, and P as functions of a for a planet 
with a lower, type B liquid water cloud that covers ~^42.3% 
and an upper, type A liquid water cloud layer that covers 
^16% of the planet for various values of the cloud optical 
thicknesses. 



are regions with only lower clouds, or only upper clouds, or 
both. 

Figure [5] shows the reflected total flux ttF^, the polar- 
ized flux Tri^n p, and the degree of polarization P for our 
model planet at A = 0.550 /im for different values of the 
cloud optical thickness, b. Because the two types of cloud 
particles have slightly different locations of the rainbows. 
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the rainbow in the total flux in the presence of two cloud 
layers appears to be somewhat broadened, as compared to 
the curves for the fluxes of single layers of clouds (the curves 
in Fig. [5] where one of the cloud optical thicknesses equals 
zero). 

In P, the strength of the rainbow feature appears to 
depend on the properties of the highest cloud layer; adding 
a lower cloud layer mainly decreases the strength of the 
Rayleigh scattering maximum by adding more low polar- 
ized light at these phase angles. P is fairly insensitive to 
the optical thickness of this lower cloud layer (i.e. 2 or 10): 
increasing the thickness of this cloud slightly increases ttFh , 
but also slightly increases 7rFn,p. This latter effect can also 
be seen in the rainbow feature in the curves for a planet 
with a single, lower cloud layer: with increasing 6, ttFh in- 
creases, but 7rFn,p increases, too! At larger phase angles 
(a > 55° in Fig.|S]), increasing h of the lower cloud layer in- 
creases Tri^n but does not increase 7rPn,p significantly, which 
can be attributed to the single scattering properties of the 
cloud droplets (see Fig. [3]). 

5. The influence of the ice cloud coverage 

On Earth, liquid water clouds can be overlaid by water ice 
clouds. These latter clouds can be entirely separate (such as 
high-altitude cirrus clouds) or the upper parts of vertically 
extended clouds that have lower parts consisting of liquid 
water droplets. The amount of ice cloud coverage on Earth 
depends strongly on the latitude. Cirrus clouds for exam- 
ple, have an annual mean coverage of ^20% in the tropics, 
while at northern mid-latitudes the annual mean coverage 
goes down to ^12%. Globally, ci rrus clouds cover ~14% o f 
the Earth's upper troposphere (jEleftheratos et al]|2007f) . 
As mentioned in Sect. (TJ water ice particles come in vari- 
ous sizes and shapes that yield single scattering properties 
that differ from those of liquid water droplets. In particular, 
light that has been singly scattered by ice cloud particles 
does not show a rainbow feature between phase angles of 
30° and 40° (i.e. between single scattering angles of 140° 
and 150°). 

To investigate the influence of water ice clouds on the 
strength of the rainbow feature of liquid water clouds, we 
used model planets with two layers of clouds: a lower layer 
of optical thickness 10, consisting of type B liquid water 
droplets and an upper l ayer contain ing water ice particles 
iHess fc Wiegneil (|l994r ): lHessl (|1998D . The lower layer is lo- 
cated between 279 and 273 K (3 and 4 km, or 0.710 and 
0.628 bar) and the upper layer between 253 and 248 K 
(7 and 8 km, or 0.426 and 0.372 bar) (c.f. Fig.[T]). The single 
scattering properties of the water ice particles clouds hav e 
been calc ulated using the code of iHess fc Wiegnerl (|1994f) : 
iHessI (|1998'). The singly scattered total and polarized flux 
and the degree of polarization of our ice cloud particles are 
shown in Fig.[3]for A ~ 0.550 /im and in Fig.|3]for A — 0.660 
and 0.865 ^m. 

Our first model planets have a liquid water cloud layer 
that covers ~ 42.3% of the surface, and an ice cloud layer 
that covers ~ 16%. About 4.5% of the liquid water clouds 
is covered by ice clouds (i.e. ~ 12% of the ice clouds is 
covering a liquid water cloud) . Figure [3] shows ttFh and P 
of our model planets as functions of a for various optical 
thicknesses of the ice clouds. 

Adding a thin ice cloud (with b = 0.5) to the cloudy 
planet lowers ttF at the smallest phase angles (cf . Fig. [5]) , 
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Fig. 9. ttFii and P as functions of a at A = 0.550 /im, for 
a planet with a lower liquid water cloud with & = 10 that 
covers ^42.3% of the planet and an upper ice cloud layer 
that covers '--^16% of the planet, and ~11% of the liquid 
water clouds. 



because the ice particles are less less strong backward scat- 
tering than the liquid water droplets (cf. Fig. [3]). At inter- 
mediate phase angles (around a = 90°), the thin ice clouds 
brighten the planet, because at those angles, their single 
scattering phase function is higher than that of the liq- 
uid water droplets (cf. Fig. [S]). Clearly, with increasing ice 
cloud optical thickness the reflected flux increases across 
the whole phase angle range (which is difficult to see at 
phase angles larger than about 150°). 

The total flux nFa does not show any evidence of a 
rainbow feature, not even for the thinnest ice cloud. The 
degree of polarization P, however, clearly shows the signa- 
ture of the rainbow, even for large optical thicknesses of the 
ice cloud particles (cf. Fig. [3]). The rainbow shows up in P 
despite the overlaying ice clouds, because the ice particles 
themselves have a very low polarization signature especially 
in the phase angle region of the rainbow (see Fig. so 
the only polarized signals arise from scattering by the liq- 
uid cloud particles and the gas molecules (which yields the 
maximum in P around 90°). Increasing the ice cloud opti- 
cal thickness, decreases P across all phase angles because it 
increases the amount of mostly unpolarized reflected light. 

In Fig. [TUl we have plotted TrPn, 7rPn,p and P for the 
same model planet as used for Fig. [HI except for an ice cloud 
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Fig. 10. Similar to Fig. [31 except for an ice cloud with 
b 2.0, and for A = 0.550, 0.660, and 0.865 ^m. The nFn 
curves for 0.660 and 0.865 /im overlap. 

optical thickness of 2.0, and for A = 0.660 and 0.865 /im. 
For comparison, the A — 0.550 fim curves (see Fig. ^ are 
also shown. With increasing wavelength, the total reflected 
flux decreases, mainly because the Rayleigh scattering opti- 
cal thickness above and between the clouds decreases. The 
decrease of the Rayleigh scattering optical thickness is also 
apparent from the decrease of the maximum in P around 
90°, and the stronger influence of the single scattering po- 
larization phase function of the liquid water cloud droplets. 



0.5 




planetary phase angle (deg) 

Fig. 11. Similar to Fig. [HI except for a lower liquid water 
cloud that covers ~25% of the planet and an upper ice cloud 
with b = 4.0. The ice cloud coverage is varied from 20% to 
75% of the water clouds. 



The ice particles leave no obvious feature in the polariza- 
tion phase function at longer wavelengths, except that they 
decrease P somewhat as compared to P of a planet with 
only a liquid cloud (cf. the 40% coverage curve in Fig. [7]). 

To investigate the influence of the ice cloud coverage 
over liquid water clouds on the reflected flux and polariza- 
tion, we used a model planet with lower liquid water clouds 
with 6 = 10 and a coverage of ~ 25%, and upper ice clouds 
with b — A right above the liquid water clouds. Figure [TT] 
shows the planet's flux and polarization phase functions at 
A = 0.550 /im for various ice cloud coverages of the liquid 
water clouds (thus, an ice cloud coverage of 75% of the liq- 
uid water clouds, equals an ice cloud coverage of ^16% of 
the planet). 

Even with ice clouds covering as little as 20% of the 
liquid water clouds (5% of the planet, but only above the 
liquid water clouds), the total flux phase function fails to 
show a hint of the rainbow. This was also apparent from 
Fig. [9j where 11% of the liquid water clouds was covered 
by an ice cloud. The polarization phase function does show 
the rainbow feature, but up from an ice cloud coverage of 
40% (a coverage of 10% of the planet) in Fig. |9l the fea- 
ture is increasingly weak such that it disappears into the 
Rayleigh scattering polarization maximum around 90°. At 
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Fig. 12. Similar to Fig. [TTl for ice clouds covering 40% 
of the liquid water clouds, and for A — 0.550, 0.660, and 
0.865 /im. 



Table 1. The cloud optical thicknesses & (at A = 0.550 /im) 
from the MODIS/Aqua data from April 25, 2011. To avoid 
having to include too many different values of 6, we have 
binned the optical thicknesses as shown below. Ice clouds 
with 6 > 20 are ignored. 
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Fig. 13. The cloud map of the Earth on April 25th, 2011, 
based on MODIS/Aqua data. The planet is covered by 
~63% liquid water clouds (gray regions) and ~36% ice 
clouds (blue regions). About 28% of the planet is covered 
by two layers of clouds. 



longer wavelengths, such as 0.865 /xm, where the Rayleigh 
scattering optical thickness above and between the clouds 
is smaller, the rainbow feature is still clearly visible in P 
when ice clouds cover 40% of the liquid water clouds. This 
can be seen in Fig. [TH 



6. Looking for the rainbow on Earth 

Figure [TT] shows us that liquid water clouds on a planet 
can be detected using polarization measurements of the 
reflected starlight, even when the liquid water clouds are 
(partly) covered by water ice clouds. The question arises 
whether the rainbow would be visible for a distant observer 
of the Earth. As far as we know, the only polarization ob- 
servations of the Earth have been done by the POLDER 
instrument (POLarization and Directionality of the Earth's 
Reflectances) a version of which is currently flying onboard 
the PARASOL satellite (for a description of the instrument, 
see iDeschamps et al.lll993 ). Locally, POLDER has indeed 
obs erved the primary r ainbow above liquid water clouds 
Csee lGoloub et al.l l2000^. However, since PARASOL is in a 
low-Earth-orbit, the POLDER measurements are not rep- 
resentative for polarization observations of the whole (disk- 
integrated) Earth observed from afar. 



In the absence of real polarization observations, we have 
simulated the total flux and polarization of light reflected 
by the whole Earth using cloud properties derived from 
observations by MODIS (Moderate Resolution Imaging 
Spectroradiometer), onboard NASA's Aqua satellite. We 
used MODIS' cloud coverage (the horizontal distribution 
of the clouds), cloud thermodynamic phase (liquid or ice), 
and the cloud optical thickness as measured on April 25th, 
2011, to build a model Earth. To limit the number of (time- 
consuming) calculations, we binned the measured optical 
thicknesses according to Table [TJ Ice clouds with an op- 
tical thicknesses larger than 20 (at A = 0.550 /zm) were 
ignored, so as to include only cirrus/cirrostratus ice clouds 
in our sample (according to the ISCCP categorization) and 
to avoid deep convection clouds. We assume that the liq- 
uid water clouds consist of the typ e B droplets and f or the 
ice particl es, we use the models of iHess fc Wieeneil (|1994[) 
and Hes^ (|1998f ). Finally, our ice clouds are positioned at 
altitudes with temperatures lower than 253 K, so that we 
avoid mixed-phase clouds (in which droplets and crystals 
exist side-by-side). The surface is assumed to be black as 
an ocean. 

Figure [13] shows the cloud map of our model Earth. On 
April 25th, 2011, about 85% of the planet was covered by 
clouds (liquid and/or ice), about 14% was covered by both 
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Fig. 14. TrFn and P as functions of a for the model Earth 
with clouds as shown in Fig. [T31 



liquid and ice clouds, about 63% of the planet was covered 
by liquid water clouds and about 36% by water ice clouds. 

In Fig. [TU we show the calculated ttF^ and P of our 
model Earth at A = 0.550 /im and 0.865 /im. The total re- 
flected flux at a = 0° equals the planet's geometric albedo, 
which is very similar at the two wavelengths and equals 
~0.22 at A = 0.55 /xm. This value is slightly smaller than 
the o ne found in l iterature for Earth's geometric albedo 
(0.33 (|Brow3[200^) and is due to our use of a black sur- 
face on our planet. It is clear that the total reflected flux 
does not show a rainbow feature, while the degree of po- 
larization does. The maximum P in the rainbow feature is 
0.08 (8% in polarization) at A=0.55 fiui, and the absolute 
difference with the nearby local minimum around a = 50°, 
is 0.02 in P (an absolute difference of 2% in polarization). 
At A = 0.865 /im, the rainbow feature is even more pro- 
nounced in P, with an absolute difference with its surround- 
ings of more than 4 %. Interestingly, P reaches zero around 
the rainbow feature (which indicates a change of direction 
of the polarization), which should facilitate measuring the 
str ength o f the r ainbow feature. 

iBailevI ()2007t ) derived the disk-integrated polarization 
signature of the partly-cloudy Earth from remote-sensing 
data and predicted degrees of polarization of ^12.7% to 
~15.5% in the rainbow peak at wavelengths between about 
0.5 to 0.8 microns. These values are higher than our val- 



ues of 4% to 8% as shown in Fig. [TH The differenc es are, 
however, small when considering the estimations bv iBailevI 
(2007) on the influence of multiple scattering and ice clouds 
on the polarization. 



7. Summary and conclusions 

The flux and degree of linear polarization of sunlight that 
is scattered by cloud particles and that is reflected back to 
space depend on the phase angle. Close to a phase angle 
of 40°, the flux and degree of polarization of unpolarized 
incident light that is singly scattered by spherical water 
cloud droplets, shows the enhancement that is known as 
the primary rainbow. The multiple scattering of light within 
clouds dilutes the primary rainbow feature both in flux and 
in polarization. In Earth-observation, with a spatial reso- 
lution of a few kilometers, the detection of the polarized 
rainbow feature is used t o discriminate betw een liquid wa- 
ter clouds and ice clouds (jGoloub et al.ll2000l ). Light that is 
scattered by ice cloud particles, such as hexagonal crystals, 
does not show the rainbow feature. Knowledge of the cloud 
properties (thermodynamical phase, optical thickness, mi- 
crophysical properties etc ) are crucial for stu dies of global 
clima te change on Earth ( Goloub et al.ll2000() . 

In lKaralidi et al.l (j201lh . we used horizontally homoge- 
neous model planets to investigate the strength of the rain- 
bow feature in flux and polarization for various model atmo- 
spheres of terrestrial exoplanets. In that paper, we briefly 
discussed the influence of a second atmospheric layer con- 
taining clouds on top of another one and the implications 
that that could have on our ability to characterize the plan- 
etary atmosphere. We concluded that even when the upper 
cloud has a relatively small optical thickness, the polariza- 
tion signal of the planet is mainly determined by the prop- 
erties of the upper cloud particles. Here, in Sects. |4] and [5] 
we extended this research on horizontally inhomogeneous 
exoplanets and for various cloud cases. 

We noticed that in case the cloud layers contain clouds 
of similar nature (for example liquid water clouds) and 
when the cloud particles present a size stratification with 
altitude it is the top cloud layer that will define the to- 
tal pla netary signal, as was also the case in iKaralidi et al.l 
(|201lD . On the other hand if the clouds are made out of the 
same nature and size particles it is the (optically) thick- 
est clouds that will define to a largest extent the planetary 
signal. 

In case the upper cloud layer contains ice clouds, we 
noticed that the characteristics of the planetary signal can 
depend on either one of the cloud layers, depending on their 
overcast and the optical thickness of the ice cloud layer. 
Using our homogeneous planet code we have seen that for 
^icc ^ 3 the existence of any lower cloud on the observed 
planetary pixel will be masked. When the overcast of the 
two cloud layers is small the existence of the ice cloud layer 
does not seem to be able to cover the existence of the water 
cloud layer (i.e. the rainbow), even for high values of the 
optical thickness (see Fig. IH]). 

Even in case the ice cloud layer has an optical thickness 
large enough to mask the existence of the underlying liquid 
water cloud, the rainbow feature of the liquid water cloud 
survives for the case the ice cloud layer covers slightly more 
than half the water clouds (see Fig. [TT]) . So unless our ob- 
served exoplanet contains a very large number of ice clouds 
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in its atmosphere, the rainbow of the water clouds will still 
be visible in the planetary P signal. 

An interesting test-case for the detection of the rain- 
bow is our own Earth. To test whether a distant observer 
would be able to detect a rainbow due to Earth's liquid 
water clouds partly covered by water ice clouds, we mod- 
eled the Earth's cloud coverage using MODIS/Aqua data 
from April 25, 2011. These data contained the location 
and optical thickness of ice and water clouds across the 
planet. We binned the cloud optical thicknesses in a limited 
number (5) of optical thickness values (to avoid too many 
time-consuming computations), as tabulated in Table [TJ 
and modelled the scattering properties of the liquid water 
cloud droplets using Mie-scatte ring and tho se of the water 
ice crystals using the models of lHessI ([l99l)- 

Our calculations for the disk-integrated flux and po- 
larization signals of this model Earth as functions of the 
planetary phase angle and at A = 0.550 fim and 0.865 fim, 
show that the flux does not have the primary rainbow fea- 
ture. Flux observations as a function of the phase angle 
would thus not provide an indication of the liquid water 
clouds on Earth. In the polarization signal, however, the 
rainbow is clearly visible, especially at the longer wave- 
lengths (A = 0.865 ^m). Polarimetry as a function of the 
planetary phase angle would thus establish the existence of 
liquid water clouds on our planet. 

The results presented in this paper were for a black plan- 
etary surface. A bright surface that reflects light with a low 
degree of polarization would not significantly affect the re- 
flected polarized flux, but it would increase the total flux, 
and hence decrease the degree of polarization. If a planet 
were thus covered by a bright surface and few clouds, the 
rainbow feature would be less strong than presented here 
(with increasing cloud coverage, the influence of the surface 
would decrease), although P at other phase angles would 
also be subdued. For example, comparing our results for a 
planet with a mean Earth cloud coverage completely cov- 
ered by a sandy surface with an albedo of 0.243 (at 0.55 
/im), polarization in the rainbow would be about 15.39%, 
compared to 41.7% for a planet completely covered by a 
black surface. Detailed calculations for planets with real- 
istically inhomogeneous surfaces and inhomogeneous cloud 
decks, preferably including the variations of the cloud deck 
in time, would help to study the effects of the surface re- 
flection. With an ocean surface with waves, one could also 
expect the glint of sta rlight to contribute to th e reflected 
flux and polarization (jWilliams fc GaidosI l2008l) . How of- 
ten this glint would be visible through broken clouds and 
its effect on the rainbow of starlight that is scattered by 
cloud particles, when it is indeed visible, will be subject to 
further studies. 

Summarizing, the primary rainbow of starlight that has 
been scattered by liquid water clouds should be observ- 
able for modest coverages (10% - 20%) of liquid clouds, 
and even when liquid water clouds are partly covered by 
water ice clouds. The total flux of the reflected starlight 
as a function of the planetary phase angle does not show 
the rainbow feature due to the presence of ice clouds. The 
degree of linear polarization of this light as a function of 
the phase angle will usually show the rainbow feature, even 
when a large fraction (up to ~50%) of the liquid water 
clouds are covered by ice clouds. Polarimetry of starlight 
that is reflected by exoplanets, thus provides a strong tool 



for the detection of liquid water clouds in the planetary 
atmospheres. 
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